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Abstract. Similar formalisms have been independently developed in 
psychology, to deal with the issue of selective influences (deciding which 
of several experimental manipulations selectively influences each of sev- 
eral, generally non-independent, response variables), and in quantum 
mechanics (QM), to deal with the EPR entanglement phenomena (de- 
ciding whether an EPR experiment allows for a "classical" account). The 
parallels between these problems are established by observing that any 
two noncommuting measurements in QM are mutually exclusive and 
can therefore be treated as analogs of different values of one and the 
same input. Both problems reduce to that of the existence of a jointly 
distributed system of random variables, one variable for every value of 
every input (in psychology) or every measurement on every particle in- 
volved (in an EPR experiment). We overview three classes of necessary 
conditions (some of them also sufficient under additional constraints) for 
the existence of such joint distributions. 

Keywords: Bell-CHSH-Fine inequalities, cosphericity test, EPR paradigm, 
joint distribution criterion, linear feasibility test, non-commuting mea- 
surements, pseudo-quasi-metrics on random variables, quantum entan- 
glement, selective influences. 



1 Introduction 

Given a set of inputs into a system and a set of stochastically non-independent 
outputs, what is the precise meaning and means of ascertaining that a given 
output is not influenced by a given input? This paper reviews the developments 
related to this question. 

The problem can be illustrated on the following diagram of selective influ- 
ences: 



a 1 = { 


M,x,y} 


a 2 = 


{x} 


a 3 = 


w,z} 
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A , A 2 , and A 3 here are random outputs, w,x,y,z are inputs (usually referred 
to as external factors in psychology and as measurement settings in QM) , and 
arrows indicate the relation "may influence": thus, the diagram does not say 
that A 2 is necessarily influenced by x, but rather that A 2 is not influenced 
by w,y,z. The diagram is shown in the canonical form, i.e., the inputs are 
redefined, {w,x,y} into a , {x} into a 2 , etc., so that each output A 1 may only 
be influenced by a single input a 1 that may not influence other outputs. We say 
then, for brevity, that (A 1 , A 2 , A 3 ) are selectively influenced by (a 1 , a 2 , a 3 ) and 
write this as 

(A 1 , A 2 , A 3 ) ^ (a 1 , a 2 , a 3 ). (2) 

Inputs (a 1 , a 2 , a 3 ) are treated as deterministic quantities, i.e., even if they are 
random variables, the joint distribution of the outputs is always conditioned 
on their specific values. Each input can have one of several values, and the joint 
distribution of {A 1 , A 2 , A 3 ) is known for each allowable treatment, a combination 
of input values. Thus, if w, x, y, z are all binary, then a 1 , a 2 , a 3 may be viewed 
as inputs with 8, 2, and 4 values, respectively, but the number of allowable 
treatments cannot exceed 16 < 8 x 2 x 4. It can be less than 16 because some of 
the combinations may be physically impossible or simply not used or observed. 

As a motivating example, consider a double-detection experiment in which 
two stimuli, say brief flashes, are presented simultaneously (right-left) or in a 
succession (first-second), each on one of two levels of intensity. The observer is 
asked to state, for each of the two observation areas (i.e., locations or time in- 
tervals), whether it contains a flash (Yes/No). The results of such an experiment 
are statistical estimates of 16 probabilities 



p(A\A 2 \a\a 2 ) = Pr 
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where a 1 (i = 1, 2) is the input representing the ith observation area (with values 
a\,a\), and A 1 is the response (Yes or No) to the ith observation area. Assume 
that A 1 and A 2 for a given (a^a 2 ) are not independent (due to attention 
fluctuations, perceptual learning, fatigue, etc.) In what sense then can we say 
that (A 1 , A 2 ) +-P (a 1 , a 2 ), and by what means can we find out if this is true? 

Many empirical situations have precisely the same formal structure. In QM, 
an example is provided by the Bohmian version of the EPR paradigm [3] : two 
subatomic particles are emitted from a common source in such a way that they 
retain highly correlated spins as they run away from each other. An experiment 
may consist, e.g., in measuring the spin of electron 1 along one of two axes, 
a\ or a>\, and (in another location but simultaneously in some inertial frame 
of reference) measuring the spin of electron 2 along one of two axes, a\ or a 2 . 



The outcome of a measurement on electron 1, A 



random variable with 



two possible values, "up" or "down," and the same holds for A 2 , outcome of a 
measurement on electron 2. The question here is: for i = 1,2, can we say that 
A 1 may only depend on a % , even though A 1 and A 2 are not independent? What 
makes this situation formally identical with the double-detection example is that 
the measurements along different axes, a\ and a\, are noncommuting, i.e., they 
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cannot be performed on the ith particle simultaneously. This makes it possible to 
consider them (measurements performed, not to be confused with their recorded 
outcomes) as mutually exclusive values of input a 1 . The results of such an ex- 
periment are described by (3), with Yes/No interpreted as spin up/down. In the 
original EPR paradigm [14] the non-commuting measurements are those of mo- 
mentum and location, each with a continuum of possible values. Our parallel with 
the issue of selective influences requires that the measurements of the momen- 
tum and of the location of a given particle be interpreted as mutually exclusive 
values of one and the same input, "(measurement of the) momentum-location of 
the particle." This may be less intuitive than the analogous interpretation of the 
spins along different axes. 

The question of selective influences cannot generally be decided based on 
the marginal distributions of the outputs alone. The most important example 
here is the classical CHSH experiment [4] where the marginal distributions of A 1 
and A 2 (in the case of two electrons) remain constant, with Pr [spin up] = 1/2. 
Examples from psychology are also readily available, especially if one adopts 
a copula view of the joint distributions. Thus, a 1 and a 2 may represent two 
stimuli presented in a succession (each having several values), and A 1 , A 2 be 
response times quantiles. The marginal distributions then are always the same, 
unit-uniform. 

2 A Historical Note 

The issue of selective influences was introduced to psychology in Sternberg's in- 
fluential paper [22] , in the context of studying consecutive "stages" of information 
processing. Sternberg acknowledged that selective influences can hold even if the 
durations of the stages are not stochastically independent, but he lacked math- 
ematical apparatus for dealing with this possibility. Townsend [24] proposed to 
formalize the notion of selectively influenced and stochastically interdependent 
random variables by the concept of " indirect nonselectiveness": the conditional 
distribution of the variable A 1 given any value a 2 of the variable A 2 , depends on 
a 1 only, and, by symmetry, the conditional distribution of A 2 at any A 1 = a 1 
depends on a 2 only. Under the name of " conditionally selective influence" this 
notion was mathematically characterized and generalized in [5]. Thus, if all com- 
binations of values of inputs a 1 , a 2 are allowable and random outputs A 1 , A 2 are 
discrete, the diagram (A 1 , A 2 ) c T- d (a 1 , a 2 ), where <— means "is conditionally 
selectively influenced," holds if and only if Pr [A 1 = a 1 , A 2 = a 2 | a*, a 2 ] can be 
presented as 

/12 (a 1 , a 2 ) h (a\al) f 2 (a 2 , a 2 ) / a 2 ) , (4) 

for all values (a 1 , a 2 ) of (A 1 , A 2 ) at all treatments (a*, a 2 ). Conditional selec- 
tivity is a useful notion, but it is not a satisfactory formalization of the intuitive 
notion of selective influences. The reason is that (A 1 , A 2 ) C T A (a 1 , a 2 ) can be 
shown [5] to violate the following obvious property of an acceptable definition: 
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the marginal distributions of A 1 and A 2 do not depend on, respectively, a 2 and 
a 1 ("marginal selectivity" [25]). 

A different approach to selective influences, reviewed below, is based on 
[6,7,9,10,11,12,19]. As it turns out 1 this approach parallels the development in 
QM of the issue of whether an EPR experiment can have a "classical" explana- 
tion (in terms of non-contextual local variables). The Joint Distribution Criterion 
which is at the heart of this development (see below) was indirectly introduced 
in the celebrated work of Bell [2], and explicitly in [15,16,23]. 

3 Basic Notions 

Aimed at providing a broad overview of concepts and results, the content of this 
paper partially overlaps with that of several previous publications, especially 
[11,12,19]. 

Random variables are understood in the broadest sense, as measurable func- 
tions X : V s — >■ V, with no restrictions on the sample spaces (V s , £ s , /it s ) and 
the induced probability spaces (distributions) (V,£,^i). In particular, any set 
X of jointly distributed random variables (functions on the same sample space) 
is a random variable, and its distribution (V, £, fx) is referred to as the joint 
distribution of its elements. We use symbol ~ in the meaning of "has the same 
distribution as." A random variable in the narrow sense is a special case of a 
random entity, with V a finite product of countable sets and intervals of reals, 
and £ the smallest sigma-algebra containing the corresponding product of power 
sets and Lebesgue sigma-algebras. Note that a vector of random variables in the 
narrow sense is a random variable in the narrow sense. 

Consider an indexed set a — {a x : X G A^ , with each a x being a set referred 
to as a (deterministic) input, with the elements of {A} x a x called input points. 
Input points therefore are pairs of the form x — (A, w), with w G a x , and should 
not be confused with input values w. A nonempty set <P C Hag^i a> ~ 1S ca U e d a set 
of (allowable) treatments. A treatment therefore is a function <$> : A — > {J\ £j \Ct x 
such that <j> (A) G a x for any A G A. 

Let there be a collection of sets of random variables A^ (A G A, <p S <P), 

referred to as (random) outputs, with distributions (V\£ x ,^y Let 

A^ = {A x : A e A} , <P e #, (5) 

be a random variable with a known distribution (the joint distribution of all A^ 
in A<j,) for every treatment We define 

A x = {A x : G , A G A, (6) 

with the understanding that A x is not generally a random variable, i.e., A^ for 
different <fi are not necessarily jointly distributed. The definition of the relation 

{A x : Xe A} ^ {a x : Xe A} , (7) 

1 This was first pointed out to us by Jerome Busemeyer (personal communication, 
November 2010), for which we remain deeply grateful. 
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interpreted as "for each A G A, A x may be influenced by a x only," can be given 
in three equivalent forms: 

(SI i ) there are independent random variables C, {S x : A G A}, and functions 

{R x (w,C,S x ) : w G a x ,X G A} , (8) 

such that, for any treatment <j) E <P, 

{R x (0(\) 7 C,S x ) : Ae/1} (9) 

(Sb) there is a random variable C and functions 

{P x (w,C) :w ea x ,Xe A} , (10) 

such that, for any treatment <j) E <P, 

{P x (0(\),C):\eA}~A^ (11) 

(J DC) there is a set of jointly distributed random variables 

H={H x :w£a x ,\eA} (12) 

(one random variable for every value of every input), such that, for any 
treatment <f> G <P, 

{h x (x) :AeA}~i f (13) 

The latter statement constitutes the Joint Distribution Criterion (J DC) for se- 
lective influences, and H is called the JDC (indexed) set. The proof of the equiva- 
lence [10] obtains essentially by the definition of a joint distribution, which seems 
to have been overlooked in the earlier derivations [15,16]. If A = {1, . . . , n} and 
all outputs A x are random variables in the narrow sense, then C in SI2 and 
C, S , . . . , S n in S 1 1 can also be chosen to be random variables in the narrow 
sense; moreover, their distribution functions can be chosen arbitrarily, provided 
they are continuous and strictly increasing on their domains, e.g., unit uniform 
[11]. 

Two important consequences of (7) are as follows: 

1. (nestedness) any subset A' of A, {A x : A G A'} {a x : A G A'}; in particu- 
lar, {A x : A € A'\ may not depend on inputs outside A' {complete marginal 
selectivity); 

2. (invariance with respect to input-value- specific transformations) for any set 
of measurable functions {F x (a) : w G a x , A € A, a € V x }, 

(B x : A G A) {a x : A G A} (14) 
where B x = [b x : G <p}, and B x = F x (x) (a$) . 



6 



E. N. Dzhafarov and J. V. Kujala 



These properties should be viewed as desiderata for any reasonable definition of 
selective influences. 

In QM, SI i corresponds to the existence of a "classical" probabilistic expla- 
nation. In psychology, statement SI i combined with auxiliary assumptions was 
used in [8] and [20] to analyze the representability of same-different pairwise 
discrimination probabilities by means of Thurstonian-type models in which two 
stimuli being compared are mapped into random entities (distributed in some 
hypothetical space of mental images) that in turn are mapped (deterministically 
or probabilistically) into a response, "same" or "different." Statement SI i was 
also used to analyze the response time distributions for parallel- serial networks 
of mental operations with selectively influenced components [13]. Note that the 
representation of the outputs A x as functions of the corresponding inputs a 1 and 
unobservable sources of randomness, A A -specific (S x ) and common (C), includes 
as special cases all conceivable generalizations and combinations of regression 
and factor analyses, with our term "input" corresponding to the traditional "re- 
gressor," and the term "source of randomness" to the factor- analytic "factor." 
This observation alone shows the potentially unlimited sphere of applicability of 
Sir. 

Statement SI2 (corresponding in QM to "classical" determinsitic explanation) 
and J DC turn out to be more convenient in dealing with certain foundational 
probabilistic issues [9] and for the construction of the working tests (necessary 
conditions) for selective influences [10,11,12,19]. The tests are discussed below. 

The following is a table of correspondences between the general terminology 
used in dealing with the issue of selective influences, and that of QM in dealing 
with EPR. 



Selective Probabilistic Causality (general) 


Quantum Entanglement Problem 


observed random output 


outcome of a given measurement 
on a given particle 


input (factor) 


set of noncommuting measurements 
on a given particle 


input value 


one of noncommuting measurements 
on a given particle 


joint distribution criterion 


joint distribution criterion 


diagram of selective influences 


"classical" explanation 


representation in the form SIi 


probabilistic "classical" explanation 


representation in the form SI2 


deterministic "classical" explanation 



4 Tests for Selective Influences 

Let H = {H* : w e a x ,X e yl} be a hypothetical JDC-set, i.e., a set satisfying 
(13) but not necessarily jointly distributed. Denoting 



{H$ w :Aei}=%^$, 



(15) 
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let % be a set of constraints imposed on possible distributions of H^. For in- 
stance, % may be the requirement that all be composed of Bernoulli vari- 
ables, or multivariate-normally distributed. 

A statement S (H ( j >1 , . . . ,H^ s ), with <pi, . . . ,<f> a G is called a test for the 
relation (7) under constraints H, if 

1. (observability) its truth value only depends on the distributions of , . . . , H t j >g ; 

2. (non- emptiness) it is not true for all possible distributions of H^, . . . 1 H^ 3 
satisfying %, 

3. (necessity) it is true if H is jointly distributed. 

If S (H ( p 1 , . . . , H c / >s ) is false for all distributions of H<p 1 , H^ s satisfying H un- 
less H is jointly distributed, the test is called a criterion for (7). In the following 
we assume that H. always includes the requirement of complete marginal selec- 
tivity: for any A' C A, the joint distribution of {^(yi')U(?>(yl-yl') : ^ *= ^ oes 
not depend on <f>(A — A'). If this condition is violated, (7) is ruled out trivially. 

4.1 Pseudo-quasi-distance tests 

A function d : H x H — > M is a pseudo-quasi-metric (p. q. -metric) on _ff if, for 
any H U H 2 ,H 3 e tf, 

(i) iJ 2 ) only depends on the joint distribution of (Hi,H 2 ), 

(n)d(H u H 2 )>0, 

(hi) d(H 1 ,H 1 ) = 0, 

(iv) d(H 1 ,H 3 )<d(H 1 ,H 2 ) + d(H 2 ,H 3 ). 

The conventional pseudometrics (also called semimetrics) obtain by adding 
the property cZ (-Hi ,H 2 ) = d (H 2 , Hi); the conventional quasimetrics are obtained 
by adding the property Pr [Hi = H 2 ] < 1 =>■ d(H\,H 2 ) > 0. A conventional 
metric is both a pseudometric and a quasimetric. 

A sequence of input points 

xi = (Xi,wi) ,...,xi = (Xi,wi) , (16) 

where Wi 6 a Ai for i = > 3, is called treatment-realizable if there are 

treatments cfi 1 , . . . , <f> 1 e <P (not necessarily pairwise distinct), such that 

{xi,xi} C (p 1 and {xi-i,Xi} C for i = 2, . . . , I. (17) 

If a JDC-set H exists, then for any p.q.-metric d on H we should have 

d(H^H^)=d(A^A^) (18) 

and 

d(H%-\,H£) =d(A%T\A%), i = 2,...,l, (19) 

whence 

l 

d(A^,A*[) <J2d(A%-\A%) . (20) 

i=2 
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This chain inequality constitutes a p. q. -metric test for selective influences. If this 
inequality is found not to hold for at least one treatment-realizable sequence of 
input points, selectivity (7) is ruled out [12]. 

It turns out that one needs to check the chain inequality only for irreducible 
treatment-realizable sequences x\, . . . ,xi, i.e., those with x\ ^ xi and with the 
property that the only subsequences {x^ , . . . , Xi k } with k > 1 that are sub- 
sets of treatments are pairs {xi, x{\ and {xi-\, Xi}, for i = 2, . . . , I. Inequality 
(20) is satisfied for all treatment-realizable sequences if and only if it holds for 
all irreducible sequences [12]. The situation is even simpler if <P = YixeA W x 
(all logically possible treatments are allowable). Then (20) is satisfied for all 
treatment-realizable sequences if and only if this inequality holds for all tetradic 
sequences of the form x,y,s,t, with x, s G {Xi} x a Xl , y,t G {A2} x a X2 7 x ^ ,s, 
y ± t, Ar ? A 2 [10]. 

Order- distances constitute a special class of p. q. -metrics, defined as follows. 
Let the distribution of H x G H be (V x , E x , fj, x ) . Let 

RC (J V Xl xV x \ (21) 

(\ 1 ,\ 2 )EAxA 

and let us write a < b for (a, b) G R. Let R be a total order (transitive, reflexive, 
and connected). We assume that for any (Ai, A 2 ) G A x A, Pr [H Xl ^ H X A is 
well-defined, i.e., { (a, b) : a G V Xl , 6 G F A2 , a ^ &} belongs to the product sigma- 
algebra over I7 Al and Z^ 2 . Then the function 

D (ff x \ , Hi\ ) = Pr [H*\ <H x l], (22) 

where -< is the strict order induced by ^, is well-defined, and it is a p. q. -metric 
on H, called order-distance [12]. 

As a simple example, consider the results of a CHSH type experiment with 
two spin axes per each of two entangled 1/2-spin particles. Enumerate the spin 
axes 1,2 for either particle, enumerate the two outcomes (up and down) of each 
measurement 1, 2 for particle 1 and 1', 2' for particle 2, and denote 

Pr [Hi =k,H* = I'] = Pr [a\^ = k, A 2 {iJ) = I'] = p kl{ij , (23) 

where k, I G {1, 2}. Define the order-distance Di by putting 1 ~ 1' -< 2 ~ 2', 
where ~ is equivalence induced by ^. We have then the chain inequality 

Pl2|12=Di(ff 1 1 ,flf) , . 

< D 1 (H 1 1 ,H 2 1 )+'D 1 (HlH^+D 1 (HlHi) = Pi 2 |n+p 2 i| 2 i+Pi 2 | 22 . ^ J 

Consider next a similar inequality for the order-distance D 2 defined by 1 ^ 2' -< 
2-1': 



Pll | 12 = T> 2 {H{,H%) 

< D 2 (Hl,Hl)+D 2 (HlH^+D 2 (HlHi) = Pim +p 22 \ 2l + Pll \ 22 . [ °> 
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By simple algebra, denoting 

Pr [Hi = k] = p k .\i., Pr [H] = I'] = p. iyj , (26) 
the conjunction of (24) and (25) can be shown to be equivalent to 

-1 < Pll|ll + Pll|21 + Pll|22 -Pll|12 — Pl-|2- -P-l\-l < 0. (27) 

One derives analogously 

-1 < Pll|12 + Pll|22 +Pll|21 -Pll|ll — Pl-|2- — P-l|-2 < 0, 

-1 < Pll|21 +Pll|ll +Pll|12 -Pll|22 < 0, (28) 

-1 < Pll|22 +Pll|12 +Pll|ll -Pll|21 — P-l|-2 < 0. 

The four double- inequalities (27)- (28) can be referred to as the Bell-CHSH- 
Fine inequalities [15,16], necessary and sufficient conditions for the CHSH type 
experiment to have a "classical" explanation. 



4.2 Cosphericity Tests 

Let the outputs all be random variables in the narrow sense. Denote, for any 
distinct Ai, A2 G A and any 6 <P with (Ai) = w\ and (A 2 ) = W2, 



= o XlX2 



(29) 



Cor [H£,H%\ =Cor 
where Cor designates correlation. Let ^1,^2,^3,^4 € ^ be any treatments with 

(30) 

Then, as shown in [19], if the components of H are jointly distributed, 



0i (Ai) = 02 (Ai) = wi; 0i (A 2 ) = 03 (A 2 ) = w 2 

04 (Ai) = 02 (Ai) = Ul'i, 04 (A 2 ) = 03 (A 2 ) = w' 2 . 



l J W\U)2 



< 



..A1A2 ^AiA2 A1A2 









(31) 



This is the cosphericity test for (7), called so because geometrically (31) de- 
scribes the possibility to place four points (w\, W2, w[, w' 2 ) on a unit sphere in 
3D Euclidean space so that the angles between the corresponding radius- vectors 
have cosines equal to the correlations. Note that an outcome of this test does not 
allow to predict the outcome of the same test applied to nonlinearly input- value- 
specifically transformed random variables. Due to (14), this creates a multitude 
of cosphericity tests for one and the same initial set of outputs A^ . 

In the all-important for behavioral sciences 2x2 factorial design (A = {1, 2}, 
each input is binary, and <P consists of all four possible treatments), the cospheric- 
ity test is a criterion for (A 1 , A 2 ) t-p (a 1 , a 2 ) if (perhaps following some input- 
value-specific transformation) the outputs are bivariate normally distributed for 
all four treatments [19]. 
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4.3 Linear Feasibility Test 

The Linear Feasibility Test (LFT) is a criterion for selective influences in all 
situations involving finite sets of inputs/outputs, A — {l,...,n}, with the ith 
input and ith output having finite sets of values, {1, . . . , fcj} and {1, . . . , m^}, 
respectively [11]. In other situations LFT can be used as a necessary condition 
because every set of possible values can be discretized. The distributions of 
= (H^ , . . . , Hj l n ) are represented by probabilities 

Pv[H} 1 =a 1 ,...,Hl=a n ] = Pr [A\ = a u . . . ,AJ = a n ] , (32) 

with 4> = (j 1 ,..., j n ) 6 $ and 

(ax,..., a n ) G {1, . . . , mi} x • • • x {1, . . . , m„} . (33) 

We consider this probability the [(ai, . . . , a n ) , (ji, • • • , j n )]th component of the 
mi • • • m n t- vector P (with t denoting the number of treatments in <P). The joint 
distribution of H in JDC, if it exists, is represented by probabilities 

Pr [H\ = h{...,Hl 1 =h 1 kl ,...,H? = K, ...,Hl= h n kn ] , (34) 

with 

(h\, kf ^)e{l mi } kl x ... x {1, ... , m„} fc " . (35) 

We consider this probability the (h\, . . . , h\ , . . . , h n , . . . , /i^)th component of 

the ("Tii)* 1 • • • (m n ) kn -vector Q. Consider now the Boolean matrix M with rows 
corresponding to components of P and columns to components of Q: let M (r, c) = 
1 if and only if 

1. row r corresponds to the [(ji, . . . , j n ) , (cti, . . . , a„)]th component of P, 

2. column c to the (h\, . . . , h\ i , . . . , h n , . . . , h kn )th component of Q, and 

3. h 1 n =01,.. .,hj n = a n . 

Clearly, the vector Q exists if and only if the system 

MQ = P, Q > (36) 

has a solution (is feasible). This is a linear programming task in the standard 
form (with a constant objective function). Let £ (P) be a Boolean function equal 
to 1 if and only if this system is feasible. C (P) is known to be computable, its 
time complexity being polynomial [18]. 

The potential of JDC to lead to LFT and provide an ultimate criterion for 
the Bohmian entanglement problem has not been utilized in quantum physics 
until relatively recently, when LFT was proposed in [26,27] and [1[. But the 
essence of the idea can be found in [21]. Given a set of numerical (experimentally 
estimated or theoretical) probabilities, computing C (P) is always preferable to 
dealing with explicit inequalities as their number becomes very large even for 
moderate-size vectors P. The classical Bell-CHSH-Fine inequalities (27)-(28) for 
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n = 2, fci = &2 = 2, mi = m2 = 2 (assuming that the marginal selectivity 
equalities hold) number just 8, but already for n — 2, ki = fc 2 = 2 with mi = 
?TJ2 = 3 (describing, e.g., an EPR experiment with two spin-1 particles, or two 
spin- 1 ^ ones and inefficient detectors), our computations yield 1080 inequalitiies 
equivalent to £ (P) = 1. For n = 3, fci = fc 2 = &3 = 2 and mi = m 2 = m 3 = 
2, corresponding to the GHZ paradigm [17] with three spin-!/2 particles, this 
number is 53792. Lists of such inequalities can be derived "mechanically" from 
the format of matrix M using well-known facet enumeration algorithms (see, e.g., 
program Irs at http://cgm.cs.mcgill.ca/~avis/C/lrs.html). Once such a system 
of inequalities S is derived, one can use it to prove necessity (or sufficiency) of 
any other system S' by showing, with the aid of a linear programming algorithm, 
that S' is redundant when added to S (respectively, S is redundant when added 
to S'). 
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